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Abstract 

Gravitational lensing of the cosmic microwave background by large-scale structure in the 
late universe is both a source of cosmological information and a potential contaminant of 
primordial gravity waves. Because lensing imprints growth of structure in the late universe 
on the CMB, measurements of CMB lensing will constrain parameters to which the CMB 
would not otherwise be sensitive, such as neutrino mass. 

In CMB polarization, gravitational lensing is the largest guaranteed source of B-mode 
(or curl-like) polarization. Future CMB polarization experiments with sufficient sensitivity 
to measure B-modes on small angular scales {i ~ 1000) can measure lensing with better 
sensitivity, and on different scales, than could be achieved by measuring CMB temperature 
alone. If the instrumental noise is sufficiently small (< 5 /xK-arcmin), the gravitational lensing 
contribution to the large-scale B-mode will be the limiting source of contamination when 
constraining a stochastic background of gravity waves in the early universe, one of the most 
exciting prospects for future CMB polarization experiments. High-sensitivity measurements 
of small-scale B-modes can reduce this contamination through a lens reconstruction technique 
that separates the lensing and primordial contributions to the B-mode on large scales. 

A fundamental design decision for a future CMB polarization experiment such as CMBpol 
is whether to have coarse angular resolution so that only the large-scale B-mode (and the large- 
scale E-mode from reionization) is measured, or high resolution to additionally measure CMB 
lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in 
polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave 
B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, 
and required control of beam systematics. 
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1 Gravitational lensing and CMB polarization 



1.1 Introduction 

Much of the progress in cosmology in the last two decades has been due to the well under- 
stood physics underlying the CMB anisotropy. The CMB promises to remain a gold mine for 
precision cosmology, and two new frontiers lie ahead. The first one is the primary purpose of 
this report, that is the polarized component that offers the prospects of detecting primordial 
gravitational waves and constraining recombination physics. Second, large scale structures 
between the last scattering surface and us alters the primary CMB anisotropy, through gravi- 
tational lensing (for a recent review of the theory see Other effects like the scattering off 
hot electrons in large scale structure (the Sunyaev-Zel'dovich effects), and through redshifting 
during the traverse of time- dependent potential fluctuations (the ISW effect) are relevant for 
temperature and will be mostly ignored here. In this section, we will show how those two 
frontiers actually merge when looking at CMB polarization at sub-degree angular scales. We 
will present how gravitational lensing of the polarized CMB constitutes a unique cosmological 
probe and the conceptual and practical challenges that arise. The large scale density fluctua- 
tions in the universe induce random deflections in the direction of the CMB photons as they 
propagate from the last scattering surface to us. The displacement angle is related to the 
projected surface density or, equivalently, the projected gravitational potential. This effect 
can be rewritten as a remapping of the primordial unlensed CMB field the following way: 

T{n) = T{n + V(p{n)) 
{Q±iU)in) = (g±zt/)(n + V0(n)) (1) 

where the deflection angle is expressed in terms of the gravitational potential as 

r^"' dz , , f D(zrcc) ~ D(z)\ 







H z) ^ ' ^ ' ' \ D Zr.,r)D Z 



where D{z) denotes the comoving distance to redshift z in the assumed flat cosmology and 
\l/(2;, x) is the zero-shear gravitational potential. In the Limber approximation, the power 
spectrum of is given by: 

Since the structures as described by the gravitational potential \1/ are not very correlated 
on large scales, the gravitational lensing effect is only relevant at small angular scales in 
the CMB. This fact has made CMB lensing observationally challenging so far. Nevertheless, 
using cross correlation between WMAP data and other tracers of large scale structures to 
increase the signal to noise, a detection of gravitational lensing in the CMB temperature 
has been achieved with marg inal significance, i.e. around 3a 0,1, fl. A direct detection in 
temperature is expected to be achieved soon with high significance thanks to on-going high 
angular resolution temperature surveys (e.g. ACT j^, SPT jsl, Planck 0]). 
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Figure 1: Left panel: Signal angular power spectrum for the E (dashed line) and B (solid lines) 
modes. The black solid dashed line corresponds to the lensing induced B modes for all 
the models considered. The light to dark red colored curves correspond to different r 
values, namely 0.43, 0.1, 0.01 and 0.001. The cosmological parameters used for this plot 
correspond to the WMAP5 ACDM+r best fit model [9|. Note that r = 0.43 corresponds 
to the 95% upper limit on r using this data-set. Obviously, for any allowed value of r, the 
lensing signal will dominate for £ > 200. Right panel: Redshift dependence of the two 
principal components {Zi and Z2 respectively) of the lensing potential angular power 
spectrum defined in Eq. (from [10|]). These curves illustrate the CMB polarization 
lensing sensitivity to moderate redshifts, i.e. up to z ~ 5. 



Promisingly, it was realized that the lensing of the CMB is more significant in polarization 
than in temperature This stems from the fact that the lensing effects on the CMB 
can be qualitatively understood as a smearing of the CMB acoustic peaks in the angular 
power spectrum. Since the CMB polarization has sharper acoustic peaks than temperature, 
the gravitational lensing effect is more significant in polarization than in temperature by 
approximately a factor of two. But the instrumental sensitivity required to detect the lensing 
effect in polarization is nevertheless higher than for temperature because of the weak degree 
of polarization of the CMB in the first place. 

However, the lensing of the polarized CMB presents several interesting features. First, as 
seen in Eq. ([1]), gravitational lensing lensing does not mix Q and U, it will nevertheless result 
in a mixing of the E and B modes because the transformation from (Q,U) to (E,B) is non-local 
(sl. EtI. 12 1 . In particular, E mode power will be transferred into B modes, generating in this 
way the largest guaranteed B-mode signal. This particular signal is totally independent from 
the existence of primordial B modes, i.e. the existence of tensor modes in the early universe 
as illustrated in the left panel of Fig. [TJ Since for realistic values of r, this B mode lensing 
signal is likely to dominate over the primordial one at sub-degree scales, it might limit our 



quest from primordial B mode 13|, [ij, ll5|] if not properly taken care of. The procedure of 



cleaning the lensing signal or "delensing" the B modes will be made explicit below. 
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Although a contaminant when trying to measure r, CMB polarization by itself contains 
unique cosmological information. Being sensitive to both the geometry of the universe and 
the growth of structure at moderate redshift (z<5) as illustrated in the right panel of Fig. [T], 
the CMB lensing breaks the angular diameter distance degeneracy in the CMB. It gives 
us a unique handle on the universe expansion history between recombination and moderate 
redshifts that is a rare probe of early dark energy. It provides access to the deepest two 
dimensional mass maps possible, thus anchoring tomographic studies of the evolution of dark 
energy at lower redshifts. CMB polarization lensing also provides a unique opportunity to 
map the distribution of matter on large scales and high redshifts where density fluctuations 
are still in the linear regime and are thus robust cosmological probes. On smaller scales, CMB 
lens reconstruction can directly probe halo mass profiles, without any need to calibrate cluster 



masses against other observables such as SZ temperature [16|, Il7|, Il8|, Il9|, l20|, l21|, |22|, |23|, [24 . 
Furthermore, since the lensing B-modes allow for an order of magnitude extension to smaller 
scales of the lensing potential as compared to temperature lensing, it is uniquely sensitive 
to p arameters that affect structure formation in the late universe, such as neutrino masses 



10, 25 



It must be said however that holding these promises is observationally demanding. Grav- 
itational lensing of the polarized CMB is a small scale manifestation of the very large scale 
properties of the intervening mass distribution. It therefore requires both high angular reso- 
lution (< 10 arcmin) and wide-field surveys (> square degrees) to be exploited. This comes of 
course at an additional cost and complexity for a satellite mission that must be quantitatively 
weighted against the scientific returns. This section aims at providing the science elements 
relevant to this debate. 



1.2 Lens reconstruction and delensing 



The most powerful techniques for extracting the gravitational lensing signal from the CMB 
are based on the idea of "lens reconstruction" , in which the defiection operation in Eq. ([1]) is 
inverted statistically: starting from the lensed (observed] CMB, one defines an estimator 



for the lens potential (which is not directly observable) [l3|, M, M, M, M, M, M, l30|, ImI, l32 , 



33|, 13J, 135|, m 



To understand intuitively how this is possible, imagine that both the lensed E-mode and 
B-mode have been measured with high signal-to-noise. Because there is no unlensed B-mode, 
the defiection operation (Eq. ([T])) converts two unobserved fields (the unlensed E-mode and 
the lens potential) into two observed fields (the lensed E-mode and B-mode). Inverting the 
defiection operation, to recover the unobserved fields from the observed ones, is possible (at 
least at the level of counting degrees of freedom) because it amounts to solving for two free 
fields given the values of two observed fieldsJl] 



^This intuitive description fails to capture some qualitative features of lens reconstruction; for example 
that lens reconstruction can be done (at lower signal-to-noise) from CMB temperature alone, or that joint 
estimation of a gradient and curl mode in the deflection angles is possible. However, it does give a simple 
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On a technical level, lens reconstruction is possible because the B-mode generated by 
gravitational lensing is highly correlated to the E-mode, with a correlation whose "shape" 
depends on the realization of the lens potential 0. In a fixed realization of the lens potential, 
the EB two-point function is of the form (see App. Rl) : 



where we have used the notation (■)cmb to emphasize that the expectation value is taken over 
CMB realizations in a fixed realization of 0. (Notation in Eq. (jl]) and elsewhere in the paper 



follows Dvorkin & Smith, to appear [37 1 



By summing (with minimum variance weighting) over two-point terms in the CMB which 
average to a given mode (p^m of the lensing potential, we can write down an estimator for 
the mode: 



ilinii2m2 



'+ 1 ^ 



IrEB 12 



This estimator is unbiased, in the sense that: 



^im / — <+>lm 

CMB 



(6) 



(7) 



and its covariance is giver|^ by: 



We therefore interpret (pim as a noisy reconstruction of the lens potential 0, with noise power 
spectrum given by the quantity iV/"^ defined in Eq. (Note that the expectation value (■) 
in Eq. ([H]) is taken over realizations of the CMB and lens potential.) 

In Fig. [21 we show some example signal and noise power spectra for the EB quadratic 
estimator, for low- noise (1 /iK-arcmin noise, 1 arcmin beam) and high- noise (10 /xK-arcmin 
noise, 10 arcmin beam) polarization measurements. As the instrumental sensitivity varies 
over this range, the lens reconstruction goes from having signal-to-noise < 1 on all angular 
scales, to being a high signal-to-noise reconstruction out to sub-degree scales (£ < 1000). For 
comparison, we also show a noise power spectrum for the TT quadratic estimator, assuming 
cosmic variance limited observations out to £max = 3000 1^ It is seen that the signal-to- noise 



intuitive interpretation of the polarization estimator in the high signal-to-noise limit. 

^Eq. ([8|) is actually an approximation; it includes most, but not all, of the contractions in the CMB four- 
point function. The additional terms can be interpreted as a change of normalization in the power spectrum 
estimator and removed using an iterative method [sHi- 

^We have chosen an ^max cutoff here, rather than assuming a noise level and beam size in temperature, 
because futuristic lens reconstruction measurements from CMB temperature are more likely to be limited by 
foregrounds on small scales than by instrumental noise (sol . 
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Figure 2: Signal power spectrum C^"^ for the CMB lens potential, and reconstruction noise power 
spectra for low- noise (1 //K-arcmin noise, 1 arcmin beam) and high-noise (10 //K-arcmin 
noise, 10 arcmin beam) polarization measurements, and for temperature measurements 
which ctrG cosmic variance limited, to ^max — 

3000. 



of the reconstruction drops sharply for i > 200 even for cosmic variance hmited observations, 
i.e. reconstructing the smallest scales in require measuring polarization and cannot be done 



from CMB temperature alone 26 



In this report, we will concentrate on two applications of the lens reconstruction estimator 
bim- First, the power spectrum estimator 



Cf" = ( ^ E $;J,Jj - Ne (9) 

V m=-l J 

is useful as a direct probe of large-scale structure, and used to constrain quantities such as 
neutrino mass to which the primary CMB is not sensitive. (Note that we define Cf*^ in Eq. 
with the noise bias term from Eq. ([8]) subtracted.) Such constraints could also be obtained 



from the B-mode power spectrum Cf^ lO, |4Q|, |41[ , but performing lens reconstruction allows 



more cosmological information to be extracted from the lensing signal. The overall signal-to- 
noise is higher, and parameter degeneracies can be broken, in cases where two parameters are 
degenerate in the B-mode power spectrum but produce distinct effects on the power spectrum 

The second application of lens reconstruction that we will study in this report is "de- 
lensing", or reducing the level of the lensing B-mode as a contaminant of the gravity wave 
signal from inflation. On an intuitive level, delensing can be described as follows. Suppose 
that the instrumental noise is sufficiently low (and foregrounds and systematics sufficiently 
well-controlled) that the lensing B-mode on large scales is the dominant source of noise when 
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estimating the tensor-to-scalar ratio (T/S). In this low-noise regime, the large-scale B-mode 
af^ has been measured with high signal-to-noise, but is a sum of lensing and primordial contri- 
butions, and sample variance of the lensing component dominates the uncertainty a(T/S). If 
we have a reconstruction (p of the lens potential, and we also have measurements of the E-mode 
on intermediate angular scales (20 ^ £ < 2000), then we can simply perform the deflection 
operation (Eq. ([1])) to obtain a reconstruction af^ of the lensed B-mode on large scales. We 
then "delens" the observed B-mode by subtracting this reconstruction {af^ —>■ af^ — af^), 
to obtain a new large-scale B-mode in which the level of lensing power has been reduced, 
while preserving the primordial contribution. An estimate of the tensor-to-scalar ratio which 
is based on this "delensed" B-mode will therefore have a smaller uncertainty a(T/S). A more 
formal version of this delensing procedure, which incorporates noise using minimum- variance 
weighting, will be given later in this report App. 



1.3 Foregrounds and systematics 

Studies of lens reconstruction to date have mainly focused on the statistical errors that can 
be obtained assuming that the microwave sky consists of a Gaussian primary CMB, lensed 
by a potential </>. In reality, there are astrophysical sources of radiation at microwave fre- 
quencies: either diffuse Galactic foregrounds (synchrotron radiation, free-free emission, dust 
emission from either vibrational or rotational modes of the grains) or extragalactic signals 
(point sources, thermal/kinetic SZ). These foreground signals are particularly worrying for 
lens reconstruction because they are not Gaussian fields, and lens reconstruction can be inter- 
preted as constraining lensing via its non-Gaussian signature in the CMB (e.g., the estimated 



power spectrum Cf' can be viewed as a trispectrum estimator [30|, l38|, |42|, |43j). At the time of 



this writing, foreground bias in lens reconstruction is largely unexplored territory (see however 



39| for some results on temperature foregrounds). In §l]we will argue that in polarization, 
the picture is relatively simple: extragalactic polarized point sources are expected to generate 
the largest foreground bias. We calculate the bias for a realistic model of the flux and redshift 
distribution of radio sources, and argue that foregrounds are not expected to bias lens recon- 
struction significantly, for a wide range of noise levels and beam sizes, if the reconstruction is 
done using polarization. 

Another practical concern for lens reconstruction from CMB polarization (or for any mea- 
surement which makes use of B-modes in a critical way) is instrumental systematics, partic- 



ulaly beam systematics [4J, |45|, |46|, |47[. Beam systematics can be classified into reducible 
(effects which are coupled to the scan strategy) and irreducible (effects which persist for an 
ideal survey), and further subclassified into specific effects (e.g. differential pointing). For 
each bea m sy stematic, the bias to lens reconstruction can be computed using the formal- 



ism from [48|, and the instrumental limit on the systematic effect (e.g. as measured from 
Jupiter maps) can be compared to the threshhold for producing a statistically significant bias 
in cosmological parameters such as rrii, or (T/S). This provides a framework for studying 
systematics that will be presented in detail in ^ 
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1.4 Outline 



The outline of this White Paper is as follows. In §21 we study CMB lensing as a source 
of cosmological information, presenting forecasts in cases where lensing adds qualitatively 
new cosmological information (compared to what could be obtained using the unlensed CMB 
alone): neutrino mass, the dark energy of state and mean curvature. In ^ we consider 
the lensing B-mode as a contaminant to the gravity wave signal on large scales, and forecast 
prospects for delensing, or reducing the level of contamination using measurements of the 
small-scale lensing potential to reconstruct the lens potential and the lensing B-mode. We 
also consider "external" delensing using datasets other than small-scale polarization: either 
small-scale CMB temperature ( §3.21) or large-scale structure ( §3.3p . but conclude that these 
approaches are not promising. In ^we consider the impact of foregrounds. We argue that 
polarized extragalactic point sources are likely to be the dominant foreground component 
for lens reconstruction, and forecast the level of contamination due using realistic modeling 
of radio sources. Finally, in ^ we study beam systematics and compute tolerance levels 
for quantities such as differential pointing or beamwidth, guided by the criterion that the 
systematic error on cosmological parameters such as (T/S) or m^, should be a small fraction 
of the statistical error. 
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2 Parameter forecasts 



In this section, we consider CMB lensing as a source of information on cosmological param- 
eters. As described in §1.21 the lens reconstruction estimator (j)£m allows us to extract a 
noisy measurement of the CMB lensing potential (pem from high-resolution observations of 
the CMB. In effect, we can observe an extra field: the resulting measurement of Cf'^ can be 
folded into a cosmological parameter analysis along with the direct measurement of the CMB 
power spectra CP, Cf^, Cf^. 



2.1 Cosmological information in the unlensed CMB 

How does the cosmological information in Cf'^ compare to the information contained in the 
CMB power spectra Cj^, Cj^ , Cf^l To answer this question, let us temporarily ignore 
CMB lensing, and ask what cosmological information is contained in the unlensed CMB 
power spectra. The qualitative picture we will give in this section is explored in much greater 
detail in e.g. 0, Hlsol, H H S . 

The shape of the CMB power spectra is directly sensitive to parameters which affect the 
physics of the evolving plasma in the early universe, such as the baryon density fift/i^, the 
matter density Vtmh'^-, and the shape of the primoridal power spectrum (parameterized through 
a spectral index Ug or perhaps additional parameters describing "running" of the spectral 
index with scale). Additionally, the overall amplitude of the power spectra is proportional to 
yl^e"^^, where Ag denotes the amplitude of the initial fluctuations and r denotes the optical 
depth to recombination. This introduces a degeneracy between and r that can be broken 
"internally" to the CMB by measuring the E-mode reionization bump on large scales, which 
is sensitive to r alone. (For more discussion of reionization and CMB polarization, we refer 
the reader to the companion white paper 541].) 

The unlensed CMB contains another parameter degeneracy, the "angular diameter dis- 
tance degeneracy", which arises when one attempts to constrain "late universe" parameters 
which mainly affect distances and growth after recombination. In this section, we will consider 
the following late universe parameters: the dark energy density f^A, dark energy equation of 
state w, curvature and neutrino mass {^rrii,). Such parameters only affect the CMB 
through the angular scale of the acoustic peaks which is a ratio of two distances: 

L = T^— (10) 

where is the angular diameter distance to recombination and s* is the sound horizon, or 
total distance that a sound wave can travel between the big bang and recombination. In a 
parameter space containing late universe parameters, one combination of the parameters 
will be well-constrained by the unlensed CMB (via the angular diameter distance -D*), leaving 
a near-perfect (A^ — l)-fold degeneracy between the others. 

Another way of describing the angular diameter distance degeneracy is that if we vary any 
of the parameters {VL^h"^ ,w,VLk}i adjusting the dark energy density I^a so that the angular 
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Figure 3: Derivatives of the unlensed power spectrum with respect to the late universe pa- 
rameters {Q.yh? ,w,Q.k} along the angular diameter distance degeneracy, showing that 
the power spectra remain constant to a good approximation. 

diameter distance to recombination D,, remains fixed, then the unlensed CMB power spectra 
remain fixed to an excellent approximation. This is directly illustrated in Fig. [3l where we 
plot the derivative of the unlensed power spectrum Cf^ with respect to each of these three 
parameters along the degeneracy surface =constant. It is seen that for cosmologically 
interesting step sizes in these parameters (say /STL^h^ = 0.01, Aw = 0.2, AQ^ = 0.01) the 
fractional change in Cf^ is very small and the unlensed CMB is essentially unchanged. 

2.2 Cosmological information from CMB lensing 

Lens reconstruction presents the possibility of breaking the angular diameter distance de- 
generacy in the unlensed CMB, by measuring the power spectrum Cf*^ of the lens potential. 
This power spectrum can be written as a line-of-sight integral which includes both geomet- 
ric distances and the power spectrum of the evolving potential (Eq. ([2])), so it depends on 
both distances and growth and is generally sensitive to late universe parameters such as 
{Qyh"^ ,w,flK}- This can be seen explicitly in Fig. |U where we show the derivative of the 
power spectrum Cf'^ with respect to each of the three parameters, taking the derivative along 
the degeneracy surface =constant as in Fig. [3J Comparing the two figures, it is seen that 
measurements of the CMB lens potential do break the angular diameter distance degeneracy. 
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Figure 4: Derivatives of Cf with respect to the same late universe parameters as in Fig. [3l showing 
a large change in the power spectrum as the parameters are varied: the angular diameter 
distance degeneracy is broken by lensing. 



allowing each of these three parameters to be constrained from the CMB alone. 

Constraining late universe parameters through lensing is a future application of CMB 
experiments which measure the small-scale modes, and for experiments which measure small- 
scale polarization in particular. As remarked in the introduction, in the limit of low noise and 
high resolution, CMB polarization experiments can ultimately reconstruct the modes of (pim 
with high signal-to-noise across a wider range of angular scales {£ < 1000) than are accessible 
using CMB temperature alone. In the next few subsections, we will present forecasts for 
parameter constraints from CMB lensing, using a Fisher matrix formalism described in detail 
in App. El 

We include unlensed temperature and polarization power spectra (TT, EE, TE) in our 
analysis and include the lensing information through the deflection angle power spectrum. 
We do not use the lensed power spectra to avoid the complication of the correlation in their 
errors between different i values and with the error in Cf'^. Using the lensed spectra and 
neglecting these correlations could lead to overly optimistic forecasts 55||. A previous study 
56| found that using lensed spectra instead of the unlensed ones (plus (f)em power spectrum) 
shrunk the expected errors on w and rrii, for their version of CMBpol by about 40% and 30% 
respectively. 

We now consider neutrino mass, dark energy and curvature in turn and forecast the 
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sensitivity of CMB alone to constrain these late-universe parameters. 



2.3 Neutrino mass 



Neutrinos are a part of the standard model of particle physics and it is now known from 
neutrino oscillation experiments that neutrinos are not massless and that the three known 
mass-eigenstates are not fully degenerate. The atmospheric |57|, l58|, l59| and solar neutrino 
experiments 60, l6l|, l62] as well as experiments with man-made neutrino beams 63|, |6J] have 
measured two mass-square differences to be close to 8 x 10~^ eV^ and 3 x 10"'^ eV^. This 
implies that there must be at least one active neutrino with a mass greater than about 0.05 
eV. Fortuitously, both CMB lensing and cosmic shear experiments can get to this level of 
sensitivity 56|, |65|]. We note that the lensing experiments are sensitive to the sum of the 
neutrino masses and it is possible that the neutrinos are highly degenerate with a sum of 
masses close to or larger than 0.15 eV. 

Limits on neutrino mass. The neutrino oscillation experiments measure the mass-squared 
differences, but not the sum of the neutrino masses. The most stringent laboratory upper 
bound on absolute neutrino mass comes from tritium beta decay end-point experiments [66| 
which limit the electron neutrino mass to < 2 eV. This could improve by an order of magnitude 
in the future with the KATRIN experiment 67|. There are other proposed experiments that 
plan to get to similar sensitivity and detection limits by searching for neutrinoless double beta 
decay [68|. A Dirac mass would elude this search, but theoretical prejudice favors and the see- 
saw mechanism requires Majorana masses. Like the CMB and galaxy shear observations, these 
future tritium end-point and neutrinoless double beta decay experiments will be extremely 
challenging. 

The current large scale structure surveys (2dFGRS, SDSS) and WMAP together already 
provide powerful constraints on neutrino mass. We know that the sum of the active neutrino 



masses is less than about 0.7 eV [9|. The sum of the active neutrino masses, niiy, is related 
to their energy density Q^h"^ ^ m^/(94eV) assuming thermally populated neutrinos. As 
mentioned earlier, at the lower end, atmospheric neutrino oscillations constrain the mass of 
at least one active neutrino to be larger than about 0.05 eV. This window from 0.05 eV to 
about 1 eV can be probed with both laboratory experiments and cosmological observations. 

A change in gives rise to many effects. First, it changes the expansion rate of the 
universe. At last scattering, this leads to a change in the sound horizon and damping length 
(of the photon-baryon fluid). The change in the sound horizon shifts the position of the peaks 
and troughs in the anisotropy spectrum while the change in the damping length (relative to the 
sound horizon) changes its amplitude. Second, the presence of a relativistic or semi-relativistic 
species has an effect on the CMB even after last scattering because it causes the gravitational 
potential to change (decay) with time. The photons traversing these potential wells red-shift 
or blue-shift, and this enhances the amplitude of the anisotropy spectrum. The above two 
effects are however degenerate with other parameters, most notably the matter density. 

There is, however, a third effect that is distinct - on small scales, the presence of a 
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Figure 5: Uncertainty a{m^) on the neutrino mass as a function of beam size and noise level for 
^max = 2000 (left panel) or ^max = 4000 (right panel) using CMB lens reconstruction, 
assuming fixed w, Qk- 



massive neutrino damps the growth of structure. The net suppression of the power spectrum 
of density fluctuations is scale dependent and the relevant length scale is the Jeans length 
for neutrinos 69|, |70|, lDJ which decreases with time as the neutrino thermal speed decreases. 
This suppression of growth is ameliorated on scales larger than the Jeans length at matter- 
radiation equality, where the neutrinos cluster like cold dark matter. Neutrinos never cluster 
on scales smaller than the Jeans length today. The net result is no effect on large scales and 
a suppression of power on small scales. This effect can be used to put constraints on the 
neutrino mass using the observed galaxy power spectrum combined with CMB observations 



72] • Eisenstein et al. 73|] predicted that the primary CMB spectrum from the Planck satellite 
can measure neutrino mass with an error of 0.26 eV. 

The alteration of the gravitational potential at late times changes the gravitational lensing 
of CMB photons as they traverse these potentials. Including the gravitational lensing effect, 
the Planck error forecast improves to about 0.15 eV 
ments capable of probing down to 0.05 eV level 25|, |5| 




74( 1 . with more ambitious experi- 
Tomographic observations of the 



galaxy shear due to gravitational lensing can a achieve similar sensitivity in m,^ [65|. The 
physics in both cases is the same: gravitational lensing. However the observations and the 
associated systematics are very different. Complementary techniques are valuable since these 
measurements will be very challenging. 

In order to forecast constraints on the neutrino mass, we marginalize over the "early 
universe" parameters {Qi,h'^,Qch'^,Q\,YY[c,T,A,ns} and ^a, but assume that the parameters 
w, Qk are fixed. (We will consider joint constraints among w, Qr} in §2.5[ ) The result 

is shown in Fig. [51 We find that a satellite mission can constrain '^^rriiy, where the sum 
is taken over neutrino species u, at the 0.03-0.12 eV level depending on the noise level Ap, 
beam width ^fwhm and maximum CMB multipole i^ax used in the lens reconstruction. 
Summary of neutrino mass. The signature of a 0.1 eV neutrino in the unlensed CMB 
anisotropy spectra is small and such small masses are only detectable through their effect 
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on lensing, which comes through their influence on the gravitational potential. Future exper- 
iments like Planck will be able to statistically detect the lensing effect and thus measure or 
put upper limits on the neutrino mass. The expected 1-a error on rrii, from Planck is 0.15 
eV, while CMBpol could get down to the 0.05 eV level to measure the neutrino mass. 



2.4 Dark energy 

Dark energy affects lensing in two distinct ways. First, the presence of dark energy implies 
faster expansion and hence a decrease in the overall growth rate. Second, dark energy can 
cluster appreciably if the equation of state is not identically 1. The second effect cannot 
be modeled unless we have a microphysical description. Two simple approaches that are 
common in the literature are to (1) model dark energy as one or more scalar field(s) with 



possibly non-canonical kinetic terms (e.g., (75|, l76|, l77|, 78, 79|) and (2) model dark energy as a 
perfect fluid with a parameterized sound speed (e.g., [80|, l8l|). We will use the first approach 
with canonical kinetic terms in the following analysis. These models are collectively called 
quintessence. The effect of the dark energy density on the growth is easy to calculate for 
small scales where the clustering is irrelevant. On larger scales, where dark energy clusters 
appreciably it is no longer possible to factor the matter density fluctuations in Fourier space 
into a part that depends on time and another that depends on wave number. This was 



investigated in detail for a constant equation of state by Ma et al. 82], who found that 

— 1/2 

dark energy clusters on scales k < = 2Vqq where dark energy has been modeled as a 
scalar field Q with effective mass Vqq which is typically not much larger than 0{Hq). The 
clustering of dark energy boosts the metric perturbations and hence lensing and thus CMB 



lensing offers a way to constrain dark energy properties SSj. The primary effect is an overall 
suppression of the growth factor except on large scales. 

The CMB lensing window function is fairly broad in redshift-space. A downside of this is 
that CMB lensing will never be competitive with SNIa observations or proposed cosmic shear 
and BAO experiments as far as measuring the equation of state of dark energy is concerned. 
However, the virtue of CMB lensing is that it is an independent alternative probe of the 
acceleration of the universe. CMBpol can measure w to a precision of 0.08-0.2 depending on 
the noise level and beam size (Fig. [6]). 

The sensitivity to a broad range of redshifts also implies that CMB lensing is a unique 
probe of dark energy (more generally clustering) at 2; > 2. Note that if wx is demonstrably 
different from -1, then dark energy must cluster on (at least) large (1000 Mpc) scales and 
then the clustering properties of dark energy, say parameterized in terms of its sound speed. 



might then be measurable (e.g., |83l|). 

The broadness of the window function also implies that the CMB is sensitive to dark energy 
properties at high redshift. For the simplest quintessence models, the contribution of dark 
energy at high redshifts is negligible. However, there is no good reason to take these models 
as more than possible examples. An important question is then that of the contribution of 
dark energy to the expansion of the universe and growth of structure in the early universe. 
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Figure 6: Uncertainty a^w) on the dark energy equation of state as a function of beam size and 
noise level for ^max = 2000 (left panel) or ^max = 4000 (right panel) using CMB lens 
reconstruction, assuming fixed Q,yh'^,^K. 



There are many motivations to consider such extensions. Among the most striking concerns 
is that of the timing coincidence: why is the vacuum energy density (or scalar field potential) 
precisely small enough to just begin dominating the energy density of the universe when the 
universe grew to its present size? In this context, models with early dark energy are arguably 



more natural [79|, |8J, |85|, |86|, |87|, |88[ than simple quintessence where dark energy emerges as 



a low-redshift phenomenon. Dark energy could also be an effect that arises on horizon scales 
such as an infra-red modification to GR (e.g., (89|). Note that the dark energy does not have 



to cause the expansion of the universe to accelerate at early times. 

We do not make predictions here for early dark energy because the predictions depend 
on the models used 85|, |90|, l91[ . If the low-redshift dark energy equation of state parameters 
are constrained by other experiments (such as SNIa), then CMB lensing should be able to 
measure the (average) high-redshift equation of state at least as well as the constant equation 
of state w considered here and thus provide an unique window into the expansion history and 
growth of structure at high redshifts. 



2.5 Curvature and joint constraints 

As a further example of a parameter constraint from CMB lensing, we consider the mean 
curvature Qk- (Historically, this was the first example of a new parameter constraint from 



breaking the angular diameter degeneracy via CMB lensing SSj.) The mean curvature is 



expected to be small in most inflationary cosmologies 92j but there are interesting inflationary 
models with = 0(10"^); this is roughly the current la upper limit from combining CMB, 
BAO, and SN datasets 93|. In Fig.[7]we show la forecasts for the uncertainty a{Qk), assuming 



that rriu and w are fixed to fiducial values. A polarization satellite can obtain a{Qk) 
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Figure 7: Uncertainty cr(Ofc) on mean curvature as a function of beam size and noise level for 
Anax = 2000 (left panel) or ^max = 4000 (right panel) using CMB lens reconstruction, 
assuming fixed nii, and w. 

(few X 10"^) from the CMB alone. 

Finally, we discuss joint constraints. In Figs. O [6l [7| we have computed uncertainties on 
each of the three "late universe" parameters {m,y, w, Qk} with the other two parameters in this 
set fixed to fiducial values. One can ask, in a parameter space in which all three late universe 
parameters are floating, can they be simultaneously constrained, or are there degeneracies? To 
quantify this, we compute the 3-by-3 correlation matrix between the late universe parameters 
in the Fisher formalism with all the early universe parameters marginalized. The result is 

1 0.34 -0.82 \ 
0.34 1 -0.63 (11) 
-0.82 -0.63 1 / 

where the ordering of rows and columns is m^,w,Qk- (This matrix was computed assuming 
Ap = 1.4 /iK-arcmin and ^fwhm = 3 arcmin. The correlations are significantly different from 
zero but not so large that we would describe this parameter space as containing a degeneracy. 
The strongest correlation is between curvature and the other parameters. This is makes 
intuitive sense given the Cf"^ derivatives shown in Fig. |H where the fractional change in Cf*^ 
with respect to curvature is approximately constant and highly correlated with the derivative 
with respect to neutrino mass and w. 

In conclusion, all-sky measurements of CMB polarization with high sensitivity and res- 
olution can qualitatively add information to the unlensed CMB: using lens reconstruction, 
the neutrino mass can be constrained to roughly (T{Y2m^) = 0.05 eV, dark energy equation 
of state to roughly a{w) = 0.15, and mean curvature to roughly a{Qk) = 2.5 x 10""^. The 
precise values will depend on the noise level and beam as shown in Figs. El El El Because the 
shape of the Cf^ power spectrum is reconstructed, and the shape dependence with respect to 
each of the three parameters w, Qk} is different (Fig. H]), there are no degeneracies in 

this parameter space although the correlations between parameters are significantly different 
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from zero (Eq. (fTTj) ). 



3 Delensing the gravity wave B-mode 



Perhaps the most exciting prospect for future generations of high-sensitivity CMB polar- 
ization experiments is constraining the tensor-to-scalar ratio (T/S), by measuring B-mode 
polarization on large scales. In inflationary models, the value of (T/S) is tied directly to the 
energy scale during inflation, so that measuring this value opens a window on the physics 



that gave rise to the initial conditions of our universe (c.f. the companion white paper |93|). 
In a real experiment, the parameter uncertainty a{T/S) will receive contributions from de- 
tector noise, foreground contamination, instrumental systematics, and contamination due to 
lensing B-modes. In this section, we will consider the last of these contributions: under what 
circumstances is gravitational lensing the limiting factor in measuring (T/S), and what are 
the prospects for reducing the lensing contamination using delensing methods? Many results 



from this section have been taken from 9J], where more details will be presented. 



It is easy to compare the contributions to a{T/S) from gravitational lensing and detector 



noise. 



The contribution from polarized foregrounds is studied in the companion white papers 
95l . l96l |: we will present some analysis of instrumental systematics in ^) If we restrict 
attention to large angular scales {i < 100), then the lensing B-mode power spectrum Cf^ 
is constant to an excellent app roximation , and the statistics of the lensed B-mode can be 



991 



treated as Gaussian [97|, |98|, |99| . Therefore lensing can be simply be thought of as an excess 
source of white noise. The RMS amplitude cTj^nsing lensing B-mode on large scales is 

~ 5 /iK-arcmin; if the instrumental noise ai^st is ^ cTj^nsing' then the lensing contribution to 
a{T/S) is smaller than the noise contribution; if (Xinst ^ '^knsing' then lensing dominates. 

For experiments with ainst ^ 5 /xK-arcmin, the only possibility for reducing the level of 
lensing contamination is to use "delensing" techniques. As described in §1.21 delensing can 
be performed whenever we have a noisy template for the CMB lens potential, and noisy 
measurements of the primary E-mode on intermediate scales. We will consider several possi- 
bilities for the template (pim'- it could either be obtained "internally" from CMB polarization 
on small angular scales ( §3.ip . or "externally" from a different dataset, either small-scale 
CMB temperature ( §3.2p or observations of large-scale structure ( §3.3p . 

In each of these cases, we will present forecasts for the parameter uncertainty a (T/S) with 
and without delensing. Our forecasting methodology is presented in detail in App. [XJ but let 
us note one key point here. The effect of delensing is to change the equivalent white noise 
level of the large-scale B-mode from the value ^i^nsing ~ 5 /xK-arcmin to some smaller value 
"^(teiensed — "^kngedfl While the prccisc valucs of a{T/ S) achievable with and without delensing 
are difficult to forecast due to considerations such as loss of modes at low I due to EB mixing 



^This statement is actually empiricial; in the forecasting methodology from App.|^ we calculate a complete 
power spectrum Cf ^ for the residual lensing B-mode B' , but for all the examples in this section, we find 
that ^ is approximately constant on large scales, so that the residual B-mode can be treated as a source 
of white noise in the forecast for a{T/S). 
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from survey boundarie^j [1Q0|, llOll . Il02l . Il03l |. the ratio is simply given by: 



a{T/ S)no delensing _ ('^Icnscd) + (^inst) 



mst/ 



:i2) 



dclonsing v'^^dclcnscd/' 

For this reason, rather than presenting forecasts for a{T/S), we will forecast the ratio in 
Eq. (fT2|) . This ratio isolates the improvement in a(T/S) due to delensing alone, independent 
of the large-scale survey geometry and mode coverage. 



3.1 Delensing using small-scale polarization 

The first approach to delensing that we will consider is to construct the template (f)£m for 
the lens potential "internally" from CMB polarization, by applying a lens reconstruction 
estimator to the small-scale E and B-modes. In polarization, the quadratic estimator 
that has been discussed previously (Eq. ([5])) can be significantly improved for low noise levels 
using an iterative, likelihood-based approach 33|, |3J]. On an intuitive level, the improvement 
arises because lensed B-mode power acts as a source of noise for the quadratic estimator, but 
the estimated lens potential can be used to "delens" and reduce the level of the lensing B-mode 
for a subsequent evaluation of the quadratic estimator, leading to an iterative estimator. One 
qualitative difference between the two estimators is that if we consider the mathematical limit 
of zero instrumental noise (neglecting real- world issues such as foregrounds and systematics) , 
the quadratic estimator will have nonzero reconstruction noise, but the iterative estimator can 
reconstruct the lens potential (p and delens the B-mode perfectly. In this idealized zero-noise 
limit there is no fundamental limit to the level of (T/S) which can be detected 15|], unlike the 



case of the quadratic estimator [13|, |32| . Our forecasting methodology for delensing includes 
the improvements from using the iterative estimator, as described in App. |X1 

In Fig. [HI we show forecasts for the improvement in a{T/S) due to delensing (i.e. the 
ratio in Eq. f[T^ ) from an "internal" lens reconstruction using small-scale CMB polarization, 
for varying noise level and beam and taking ^max = 4000 throughout. Let us note some 
qualitative features of this figure. For large beam size, delensing is not very effective since it 
depends on being able to reconstruct the lens potential indirectly through its effect on the 
small-scale modes of the CMB. The effective noise level for low-£ B-modes will simply be the 
sum of lensing {cr^^^^^^ ~ 5 /iK-arcmin) and instrumental contributions. As the beam size 
decreases, delensing can improve the lensing contribution, but the instrumental contribution 
is unchanged. In the limit of a very small beam, the delensed noise level cr^j^j^^g^ will reach an 
intermediate value which is less than cxj^nsod but somewhat larger than the instrumental noise 
(T^sf For example, with 1 /iK-arcmin instrumental noise and a 2' beam, we find cr^j^j^^g^ = 1.7 
/^K-arcmin, resulting in a factor 7 improvement in a{T/S) relative to the no-delensing case, 
as shown in Fig. [HI 



""The most critical issue when forecasting a{T/S) is whether the gravity wave B-mode can be constrainted 
through the reionization bump at « 8, or whether only the recombination bump at £ 60 is measurable 
in the presence of foregrounds and sky cuts. At the level of a naive mode-counting forecast, the reionization 
bump has « 10 times the signal-to-noise of the recombination bump when constraining (T/S). 



21 



14 
12 

I 10 

5 8 

1 

5 6 

g 

g 



"O 5 10 15 20 

*Fnmi (arcmin) 

Figure 8: Ratio of a{T/S) with and without polarization delensing, forecasted using Eq. (jl2p for 
varying instrumental noise level and beam. 

Therefore, for B-mode experiments which are lensing-hmited {cr^^^ < 5 /iK-arcmin), the 
large and small scales are intimately linked: measuring the gravity wave B-mode on large 
scales ultimately depends on reconstructing the lens potential via the lensing B-mode on 
large scales. The effective noise level for constraining (T/S) has a nontrival dependence on 
the beam size as shown in Fig. [HI if the noise level and beam size are small, then large 
improvements in a{T/S) are possible. In practice, since beam size is a primary driver of cost 
and complexity (particularly for a satellite mission), the improvement in a{T/S) which we 
have forecasted in Fig. [H] will be one factor to be weighed against others when designing an 
experiment. 

3.2 Delensing using small-scale temperature 

In this subsection and the next, we consider situations in which lensing-limited (< 5 /iK- 
arcmin) CMB polarization measurements have been made on large angular scales, using an 
instrumental beam which is too large to observe the small-scale lensing B-mode needed for 
"internal" delensing. Is it possible to delens the large-scale B-mode using "external" mea- 
surements from other datasets? In order to perform delensing, we must have: 

1. A (noisy) template (pim for the CMB lens potential 

2. A (noisy) measurement Eim of the CMB E-mode on intermediate {i < 2000) angular 
scales (e.g. from the Planck satellite [3]). 

The improvement in a{T/S) which can be achieved using delensing will depend on the noise 
levels in both (pi^n and Eim- 
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Figure 9: Reconstruction noise power spectra Nf"^ from CMB temperature alone, assuming cosmic 
variance limited measurements for varying values of 



^max- 



One possible way to get the template would be to apply the quadratic estimator 
to small-scale CMB temperature measurements from another experiment with high angular 
resolution. Experiments are already underway (e.g. ACT or SPT j^) with sufficient 
sensitivity to measure the CMB temperature with high signal-to-noise, far into the damping 
tail of the temperature power spectrum (£ > 2000). For such experiments, the limiting factor 
in lens reconstruction is likely to be the presence of non-Gaussian secondary anisotropies 
(which become increasingly important as £ increases), rather than instrumental sensitivity 
or resolution. However, at the time of this writing it is unclear what ran ge o f scales will be 



"sufficiently Gaussian" to use for lens reconstruction in a real experiment [39 . 

We will model this unclear situation in an approximate way by introducing a cutoff mul- 
tipole ^max5 ^"^^ assuming that temperature multipoles I < i-^^^ can be used for lens recon- 
struction with the full statistical power of a Gaussian field (i.e. without introducing extra 
systematic error from secondary anisotropies), but multipoles i > f^g^^ are not useful for 
lens reconstruction. In Fig. [HI we show noise power spectra Nf'^ obtained using lens recon- 
struction from CMB temperature, for varying £^a^. As increases, a high signal-to-noise 
reconstruction is obtained on large scales, but on angular scales which are smaller than the 
CMB acoustic peak scale {i > 200), the reconstruction always has poor signal-to-noise. 

We would now like to forecast the improvement in a(T/S) due to delensing, i.e. the ratio 
in Eq. fll2p . In addition to the power spectrum Nf^ of the noise in the lensing template, 
this ratio will depend on the noise Nf^ in the E-mode measurement on intermediate scales 
(i.e. item 7^2 in the list above) and the instrumental noise a^^^^ on the large-scale B-mode. 
However, an upper bound on the ratio can be obtained by neglecting these noise sources and 
assuming Nf^ = a^^^^ = (but keeping the nonzero Nf'^ shown in Fig. [H]). This upper bound 
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Figure 10: Improvement in a{T/S) due to delensing (Eq. (jl2p ) from CMB temperature alone, 
assuming cosmic variance limited measurements for varying values of ^max (from [9j] ) . 



T 

max' 



is shown in Fig. [TO] for varying (f^^^ (taken from 9J]). It is seen that, even for large 
the improvement is modest: a factor of two at =3500. In practice, this will be further 
degraded by the noise sources that have been neglected in obtaining this upper bound. We 
interpret this as a negative result: it is not possible to delens CMB polarization using lens 
reconstruction from CMB temperature alone, because only the large-scale modes in can be 
reconstructed with high signal-to-noise (Fig. [9]). 



3.3 Delensing using large-scale structure 

We next consider another case of "external" delensing: using large-scale structure between 
the observer and recombination to obtain an external template (p. One could imagine using 
different f lavo r s of l arge- scale structure data (for example, cosmic shear 1041 ] or 21-cm tem- 
perature lOSl . Il06l . IIOTI ]) to construct </>, weighted to minimize the power spectrum of the 
residual field {(j) — (p), where (j) is the true CMB lens potential. 

In the previous subsection, we obtained an upper bound on the improvement in a{T/S) 



that could be obtained using temperature multipoles i < 



by making some idealizing 



assumptions: we neglected noise in the in the E-mode on intermediate scales, and in the 
B-mode on large angular scales. Since our final result showed only a modest improvement 
(Fig. (ITO!) ) even under these assumptions, we could intepret it as a general "no-go" theorem: 
lens reconstruction from CMB temperature is of very limited utility in delensing the large-scale 
B-mode. In this subsection, we will construct an analogous upper bound for the improvement 
in a{T/S) that can be obtained using measurements of large-scale structure from redshifts 

Z ■S^max* 
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Figure 11: Improvement in a{T/S) assuming perfect delensing of all large-scale structure from 
redshifts < ^max- 



If we write the CMB lens potential as a line-of-sight integral (Eq. ([2])) with contributions 
from different redshifts, then by causality alone, contributions from redshifts > Zmax cannot 
be reconstructed using large-scale structure, and must therefore be treated as "noise" power 
in the reconstruction. (The reconstruction noise is defined to be the difference between the 
true CMB lens potential and the template constructed from large-scale structure.) More 
formally, a lower bound on the noise power spectrum Nf'^ can be obtained by simply cutting 
off the redshift integral for Cf'^ at redshift ^max^ 



Hiz 



D{z,,,)D{z) 



(13) 



We now forecast the improvement in a{T/S) from delensing, making idealizing assumptions: 
we assume that the reconstruction noise is equal to the lower bound in Eq. f[T51) . and that 



cr 



B 



as in the previous subsection. The result is shown in Fig. [TT] (taken from 
94|). Since we find only a modest improvement (a factor of 2.2 for 2max =3) even with our 



idealizations, we interpret this as a general negative result: it is not possible to delens CMB 
polarization using large-scale structure. 

This conclusion assumes that only large-scale structure from redshifts ^ 3 is available 
with sufficient statistical power to construct the template (pem- One possible exception to 



this assumption may be a fu turistic 21-cm experiment such as SKA or FFTT 108| . This 



possibility is studied in |105l ] but is unlikely to be available for the next few generations of 
CMB polarization experiments. 
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4 Polarized foregrounds 



Lens reconstruction relies on the non-Gaussian nature of the lensed CMB: each mode (pim of 
the lens potential induces a small deviation from Gaussian statistics, and this permits the 
potential to be reconstructed. Because astrophysical foregrounds are non-Gaussian they are 
particularly worrisome for lens reconstruction. For example, it is not clear a priori how to 
relate the "strength" of a foreground contaminant at the power spectrum level to the bias that 



it produces in the reconstructed potential. Detailed arguments which will be presented in |109 
will show that polarized extragalactic point sources are expected to be the largest contaminant 
of the reconstructed potential. For this reason we will not analyze, e.g. polarized synchrotron 
or dust emission in this section, but we will do a detailed analysis of the contamination from 
polarized point sources. The extragalactic point sources can be modeled sufficiently well that 
a reliable estimate of the lens reconstruction bias can be made, at least at a rough level. This 



section is an abridged version of 1091], where more details will be given. 



4.1 Polarized point sources: forecasting machinery 

Naively, the contribution of point sources to reconstruction of the lensing power Cf'^ is a 
four-point function in the locations and polarization angles of the sources. However, this 
picture simplifies if the polarization angles of distinct sources are assumed to be uncorrelated: 
the only terms which generate a nonzero expectation value are 1-source and 2-source terms. 
In the absence of any observational evidence to the contrary, we will make this assumption. 




It can then be shown 109l | that the contribution of point sources to lensing reconstruction is 
contained in two effective power spectra Cf^, Cf^ which are defined by: 

S^S^YU^^iYYem'in^)^ = Cfhe^mm' (14) 

where (•) denotes an average over realizations of the point source model, summation indices 
i,i run over point sources in a given realization, and (rtj. Si) denote the location and polarized 
fiux of source i. 

The Cf power spectrum is due to the correlation of point sources with the lensing potential 
0. Extragalactic point sources are biased tracers of the large-scale matter distribution, and 
thus correlated with 0(n) through the line-of-sight lensing integral of Eq. ([2]). This leads to 
a bias in the reconstructed power spectrum given by 



AC, - -1-1 2^ A,^,^, I _ 1 (16) 
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where we have defined: 



The C^^ power spectrum, on the other hand, encapsulates the bias due to auto-correlations 
amo ng; the point sources. The relationship between the Cj^ and lensing bias is quite involved 



109|. For the purposes of this report, however, we will be neglecting the auto-clustering of 
point sources and may set {SfSj) = {Sf){Sj). In this case, C^^ has the form 

Cf = /i + /2-5.o. (18) 

We will refer to the biases originating from Ji and I2 respectively as the 1-pt and 2-pt Poisson 
terms. This form of C^^ greatly simplifies the the calculation of the bias due to Cf^, and we 
find an additive effect on lens reconstruction given by 



K<i>* A0 



1, (2^ + 1) 



2 



2 -20M2 -20/ V2 2 -4/V2 2 -4 



16vr^(2£+l)(Ci^^ + iV,f^m^^ + iV,^^)2 ^ ^ 

Our picture of the point-source bias to lensing reconstruction thus consists of two terms: 

• A "multiplicative" bias C*^*^^" due to the cross-correlation of point sources with large- 
scale-structure. 

• An "additive" bias Q'f^^'^^'' due to the non-zero four-point function of the point sources. 
This in turn separates into 1-pt and 2-pt Poisson contributions. 

Our analysis of point source bias now has two remaining steps. First, we must estimate 
the power spectra and C^^ from current observational constraints on polarized sources; 
second, we forecast the point source bias given values of these power spectra. These steps are 
carried out in the following two subsections. 



4.2 Polarized point sources: modeling 

We will estimate the point source bias at a fiducial frequency v = 100 GHz, under the following 
simplifying assumptions: 

1. We will only consider radio point source contamination, assuming that the contribution 
from infrared sources is smaller or comparable in order of magnitude at 100 GHz. 

2. As previously discussed, we will ignore point source clustering, i.e. only consider the 
one-halo contribution to C?^. 
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We expect the resulting estimate of the bias to be corre ct at the order-of-magnitude level, 
deferring a more detailed forecasting to future work 109 . 

W e wi ll model the source counts of radio sources at 100 GHz using the fitti ng fu nction 
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from [110|, based on extrapolating multifrequency observations from 15-43 GHz 

dN _ No 

Is 

where (3 = 2.15 and A'^o = 12 Jy^'^^ sr~^. Fitting functions based on independent datasets 
have been also proposed in 



(20) 
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113| and agree within a factor 2. 



We will model the polarization fraction by assuming that the Q, U components of the 
polarization are Gaussian distributed and that the polarization angle is uniform distributed. 
Under these assumptions, the PDF for the polarization is: 



dQ dU 
— ^exp 



(21) 



where we take RMS polarization fraction 7 = 0.1. This choice is somewhat conservative; 
bright sources at 20 GHz are typically 1-5% polarized but th ere i s some observational evidence 
for an increasing polarization fraction with decreasing flux 1141 . 

To calculate the flux integrals of Eqs. (IT^ [T5l) . we will take a simplified view of the point 
source removal process. We suppose that all sources above some limiting flux S^^^ have 
been identified and masked with 100% completeness, and those below are untouched. We 
further assume that this masking is performed in temperature, as it is here that our current 
understanding of CMB source extraction is best developed. Neglecting the use of polarization 
data to mask sources is an additional conservative choice in our analysis. Although the signal 
strength is weaker in polarization by a factor of 7^, the confounding CMB signal is smaller 
as well, and improvements in detection may be made. For a temperature threshold S, 
effective point source power spectra are given by 



the 
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P^s{r,k = i/D{z)) (23) 



where b denotes the (assumed constant) bias of the radio point sources, and the remaining 
notation follows Eq. Q. 

To compute the cross spectrum C^'^, we need a model for {d'^n/dS dz), the joint flux- 
redshift distribution of the point sources. We use the fitting function proposed in 115 | 
("RLF-1", fiat-spectrum) and extrapolate to 100 GHz by adding a (weakly 2;-dependent) 
normalization so that the source counts (dn/dS) are consistent with the power law in Eq. (l20l) . 



4.3 Polarized point sources: forecasts 

We may now proceed to evaluate the magnitude of point source contamination for EB lensing 
reconstruction. 
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Figure 12: Comparison of ttie signal power spectrum and the bias terms in Eqs. (j22|) . (j23p . 
for an experiment with Ap = 4 /xK-arcmin, ^fwhm = 7 arcmin, S'™'*^ = 200 mJy. 

In Fig. [121 we show the predicted bias for a (7' FWHM, Ap = 4yuK — arcmin) experiment, 
with S^°f_= 200mJy. This is the 100% completeness hmit expected for the Planck lOOGHz 



channel |116l |. and should therefore be readily achievable. For illustration, we plot the biases 
due to the components of Eqs. (l22l) . (l23l) . Their behaviour here is characteristic of all of 
our forecasts. The 1-pt poisson term has a contribution which increases rapidly with i, and 
typically dominates hj £ = 1000. The 2-pt poisson term is always subdominant, except at 
extremely low-i. The cross term C|"^ gives the largest contribution on intermediate scales. 
Note that these terms have different scalings with 5*™^^. The two poisson-type biases are 
due chiefly to point sources immediately below the removal threshold, whereas the Cf^ term 
receives contributions from all of the unresolved sources, and thus scales more slowly with 

To consider the effect of point source contamination more thoroughly we will need an 
estimate of the 5™^^ which is achievable for a given experiment. The subject of point source 



extraction is an active one, with many techniques in active development |117l ]. For summary 
purposes, however, we will take the following simplified model of this process. Suppose that 
we clean for point sources internally, by identifying all of the peaks in the CMB which are 
greater than 5a, relative to the total variance of the map. In temperature, the Ax^ due to a 
single point source with a flux of St is given by 

= cJT + Nf^- ^^^^ 
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Solving for Ax^ = 25 then gives 

( I 2£ + 1 \ 

^^'^ = ( 100^ ^ cfTivF j ■ ^^^^ 

This simple model gives values for the residual poi nt source flux in reasonable agreement with 



those determined in more complete analyses llq . In line with our other assumptions, it is 



also a somewhat pessimistic estimate. It neglects, for example, our ability to increase the 
contrast of point sources by differencing CMB maps at multiple frequencies. 

We also consider confusion as a lower limit on our ability to mask sources. Taking the 
differential number counts of Eq. (120|) . the number of point sources above a cutoff is 
given by 

]\T ( Qmax\l-/3 

N{s > srn = J_ — (26) 

To avoid confusion due to overlapping sources, we must ensure that the typical spacing be- 
tween sources is M beamwidths, for some reasonable value of M. We therefore require that 

^0 r,. n i2^^ 



^max > ^max^^^) ^ . e^^^^Yj (27) 

In what follows, we will take M > 10. Satisfying Eq. (|27l) then ensures that < (9(5%) of the 
sky will be excised for point source removal. This has the benefit of limiting the issues with 
E/B mixing due to the masking process which we have otherwise neglected. 

We would also like to establish a connection between biases at the power spectrum level 
which we calculate here and parameter constraints. For this, we use m,y as a canary parameter, 
and plot the quantity 

act 



dra^ 



cr{m,) (2^ 



If the bias to the power spectrum is well below this level, we expect that the effect on 
parameters will be negligible. 

In Fig. [13] we have plotted expected bias levels for 30', 15', and 2' FWHM experimental 
configurations. For all three beam sizes, we find that for reasonable values of S*™*^^, the 
point source bias generally small. Where it grows large enough to potentially bias parameter 
determinations, the (S/N) of the lens reconstruction is always < 0.1, and so we expect that 
the potentially contaminated high-£ reconstruction may be disposed of without significant 
loss of information. 

There are two caveats with these findings, however: 

• For a blunt (30' FWHM) beam experiment, removal of point sources to the expected 
level may be hampered by the large beam size. 

• For a sharp (2' FWHM) beam experiment, our findings are particularly dependent on 
the source count relation of Eq. fl20l) holding true to at fluxes of < 20mJy, where 
complete measurements have not yet been made. 
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Figure 13: Lensing reconstruction biases. Black/Green solid curves are CY 1^1 divided by 100. 

The black dashed curves are Am^{Cf'^) (Eq. [27]). The curves labelled in mJy are the 
total biases for the corresponding value of S^^^. Solid blue curves are for a 'fiducial' 
^max_ smaller of 200m Jy and the value determined from Eq. ([25]). The upper and 
lower dashed blue curves correspond to S^'^ for M = 20 and M = 10 (only M = 10 is 
shown for the 2' experiment for clarity, as the M = 20 curve overlaps with the fiducial 
curve). The S*™*^^ = 19.1mJy curve for the 2' experiment is the value of Eq. (|25p 
assuming that only multipoles £ < 3000 are used for source cleaning. 
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Exempting these two possible issues, to the best of our current modehng abihty it appears 
that radio point sources will not constitute a significant difficulty for a future polarization 
based lensing reconstruction. 
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5 Beam systematics 



To calculate the effect of beam systematics we invoke the Fisher information- matrix formalism. 
Our objective is to determine the susceptibility of certain cosmological parameters to beam 
systematics. We represent the extra noise due to beam systematics by analytic approximations 



47[ and include lensing extraction in the parameter inference process, following |25l . |56| for 
neutrino mass (and other cosmological parameters) inference from CMB data. 

Our main concern is the effect on the tensor-to-scalar ratio r and the total neutrino mass 
My (assuming three degenerate species). The lensing- induced B-mode signal is sensitive to 
neutrino masses and therefore a large enough beam systematic which leaks temperature or 
&mode polarization to B-mode polarization can bias the inferred neutrino mass. 

5.1 General 

Beam systematics due to optical imperfections depend on both the underlying sky, the prop- 
erties of the polarimeter and on the scanning strategy. An instructive example is the effect 
of differential pointing. This effect depends on the temperature gradient to first order. CMB 
temperature gradients at the 1°, 30', 10', 5' and 1' scales are ~ 1.4, 1.5, 3.5, 2.5 and 0.2 
/i/T/arcmin, respectively. Therefore, any temperature difference measured with a dual-beam 
experiment with a ^ 1' pointing error and non-ideal scanning strategy which is dominated 



by its dipole and octupole moments [47| will result in a ~ l/iK systematic polarization which 
has the potential to contaminate the B-mode signal. 

Similarly, the systematic induced by differential ellipticity results from the variation of 
the underlying temperature anisotropy along the two polarization-sensitive directions which, 
in general, differ in scale depending on the mean beamwidth, degree of ellipticity and the 
tilt of the polarization-sensitive direction with respect to the ellipse's principal axes. For 
example, the temperature difference measured along the major and minor axes of a 1° beam 
with a 2% ellipticity scales as the second gradient of the underlying temperature which on 
this scale is ~ 0.2/ii^/arcmin^ and the associated induced polarization is therefore expected 
to be on the ~ level. The spurious signals due to pointing error, differential beamwidth 
and beam ellipticity all peak at angular scales comparable to the beam size (since they are 
associated with features in the temperature anisotropy on sub-beam scales on the one hand 
but suffer beam dilution on yet smaller scales). If the beam size is ~ 1° the beam systematics 
mainly affect the deduced tensor-to-scalar ratio, r. If the polarimeter's beamwidth is a few 
arcminutes the associated systematics will impact the measured neutrino mass m,^, spatial 
curvature fifc, running of the scalar spectral index a and the dark energy equation of state w 
(which strongly affects the lensing-induced B-mode signal). It can certainly be the case that 
other cosmological parameters will be affected as well. 

Two other spurious polarization signals we explore are due to differential gain and differen- 
tial rotation; these effects are associated with different beam 'normalizations' and orientation, 
respectively, and are independent of the coupling between beam substructure and the under- 
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lying temperature perturbations. In particular, they have the same scale dependence as the 
primordial temperature anisotropy and polarization power spectra, respectively, and their 
peak impact will be on scales associated with the CMB's temperature anisotropy (~ 1°) and 
polarization 10'). 

5.2 The effect of systematics on lens reconstruction 

Gravitational lensing of the CMB is both a nuisance and a valuable cosmological tool (e.g. 
jsl). It certainly has the potential to complicate CMB data analysis due to the non-gaussianity 
it induces. However, it is also a unique probe of the growth of structure in the linear, and 



mildly non-linear, regimes (redshift of a few). In [25|, |56| as well as elsewhere it was shown 
that with a nearly ideal CMB experiment, neutrino mass limits can be improved by a factor 
of approximately four by including lensing extraction in the data analysis using CMB data 
alone. 

This lensing extraction process is not perfect; residual noise will affect any experiment, 
even an ideal one. This noise will, in principle, propagate to the inferred cosmological pa- 
rameters if the latter significantly depend on lensing extraction, e.g. neutrino mass, a and 
w. It is important to illustrate first the effect of beam systematics on lensing reconstruction. 
By optimally filtering the t empe rature and polarization, the lens potential can be recovered 



using quadratic estimators |118l |. It was shown that for experiments with ten times higher 
sensitivity than Planck, the EB estimator yields the tightest limits on the lens potential. This 
conclusion assumes no beam systematics which might significantly contaminate the observed 
B-mode. 

We illustrate the effect of differential gain, beamwidth, beam rotation, ellipticity and differ- 



ential pointing (see 47|) on the noise of lensing reconstruction with 20 different combinations 
of noise levels and resolution for CMBPOL spanning the sensitivity and resolution ranges 1-6 
yLiK-arcmin and 3-30 arcmin, respectively. These are perhaps the most pernicious systemat- 
ics. Beam rotation induces cross-polarization which leaks the much larger E-mode to B-mode 
polarization and differential ellipticity leaks T to B. The modified noise in reconstructing the 
lens potential, A'^;'''^, is consistently substituted into our Fisher-matrix analysis. 

5.3 Error forecast 

Accounting for beam systematics in both Stokes parameters and lensing power spectra is 
straightforward. In addition, the detector noise associated with the main beam is accounted 
for, as is conventional, by adding an exponential noise term. Assuming gaussian white noise 

« - 5^ (29) 

where a runs over the experiment's frequency bands. The noise in channel a is (assuming a 
gaussian beam) 
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where is the noise per beam in /xK, 6a is the beam width, and we assume noise from 
different channels is uncorrelated. The power spectrum then becomes 



N, 



X 



(31) 



where X is either the auto-correlations TT, EE and BB or the cross- correlations TE, TB 
and EB (the latter two power spectra vanish in the standard model but not in the presence 
of beam systematics and exotic parity- violating physics (e.g. 0, H, E3) or primordial 
magnetic fields (e.g. 121| ). For the cross-correlations, the vanish as there is no correlation 
between the instrumental noise of the temperature and polarization (in the absence of beam 
systematics) . 



5.4 Results 

We consider the effect of both irreducible and reducible systematics. By 'reducible' we refer 
to systematics which depend on the coupling of an imperfect scanning strategy to the beam 
mismatch parameters. These can, in principle, be removed or reduced during data analysis. 
This includes the differential gain, differential beamwidth and first order pointing error beam 
systematics. By 'irreducible' we refer to those systematics that depend only on the beam 
mismatch parameters (to leading order). For instance, the differential ellipticity and second 
order pointing error persist even if the scanning strategy is ideal. For reducible systematics 
the scanning strategy is a free parameter in our analysis (under the assumption it is non- 
ideal, yet uniform, over the map) and we set limits on the product of the scanning strategy 
(encapsulated by the /i and /2 parameters) and the differential gain, beamwidth and pointing. 



as will be described below. The exact definitions of /i and /2 are given in [47|, |48| but here 
we give approximate expressions under the assumption that the scanning strategy does not 
contain significant hexadecapole moment. 

h = k\B\' + \C\') (32) 



2 



where 



A ^ {exp{2ia)) 
B ^ (exp(3za)) 

C ^ (exp(m)) (33) 

and therefore while /i captures the quadrupole moment of the scanning strategy, /2 is a 
measure of its dipole and octupole moments. To calculate the power spectra we assume 
the concordance cosmological model throughout; the baryon, cold dark mater, and neutrino 
physical energy densities in critical density units Qbh"^ = 0.021, Q^h"^ = 0.111, Qi,h'^ = 0.006. 
The latter is equivalent to a total neutrino mass M^, = Yl'i=i ^v,i =0.56eV, slightly lower than 
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the current limit set b y a l oint analysis of the WMAP data and a variety of other cosmological 
probes (0.66eV, e.g. |122| ). We assume degenerate neutrino masses, i.e. all neutrinos have 
the same mass, 0.19 eV, for the purpose of illustration, and we do not attempt to address here 
the question of what tolerance levels are required to determine the neutrino hierarchy. As 
was shown in 25|] , the prospects for determining the neutrino hierarchy from the CMB alone, 
even in the absence of systematics, are not very promising. This conclusion may change when 
other probes, e.g. Ly-a forest, are added to the analysis. Dark energy makes up the rest of 
the energy required for closure density. 

We limit our analysis to the tensor-to-scalar ratio r, and total neutrino mass M,^. While 
r is mainly constrained by the primordial B-mode signal that peaks on degree scales (and 
is therefore not expected to be overwhelmed by the beam systematics which peak at sub- 
beam scales), it is still susceptible to the tail of these systematics, extending all the way 
to degree scales, because of its expected small amplitude (less than O.lfiK). The tensor-to- 
scalar ratio is also affected by differential gain and rotation which are simply rescalings of 
temperature anisotropy and E-mode polarization power spectra, respectively, and therefore 
do not necessarily peak at scales beyond the primordial signal. 

Ideally, the lensing signal, which peaks at / ~ 1000, provides a useful handle on the neu- 
trino mass as well as other cosmological parameters which govern the evolution of the large 
scale structure and gravitational potentials. However, the inherent noise in the lensing recon- 
struction process 26| which depends, among other things, on the instrument specifications 
(detector noise and beamwidth), now depends on beam systematics as well. The systematics, 
however, depend on the cosmological parameters through temperature leakage to polarization, 
and as a result there is a complicated interplay between these signals and the information 
they provide on cosmological parameters. As our numerical calculations show, the effect on 
the inferred cosmological parameters stems from both the direct effect of the systematics on 
the parameters and the indirect effect on the noise in the lensing reconstruction, Nf'^, in cases 
where the MV estimator is dominated by the EE correlations. 

The Fisher information-matrix gives a first order approximation to the lower bounds on 
errors inferred for these parameters. We follow 46| in quantifying the required tolerance on 
the differential gain, differential beamwidth, pointing, ellipticity and rotation. To estimate the 
effect of systematics and to set the systematics to a given tolerance limit one has to compare 
the systematics-free la error in the i-th parameter to the error obtained in the presence of 
systematics. The latter has two components; the bias and the uncertainty (which depends 
on the curvature of the likelihood function, i.e. to what extent does the information matrix 
constrain the cosmological model in question). As in 46| we define 

I A? 



/3 = ^L- (34) 

^A, ' 

where the superscript refers to values evaluated at the peak of the likelihood function, 
i.e. the values we assume for the underlying model, and AAj and Acta^ are the bias and the 
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change in the statistical error for a given experiment and for the parameters Aj induced by 
the beam systematics, respectively. As shown in 46|] these two parameters depend solely on 
the primordial, lensing and systematics power spectra. We require both 6 and /3 not to exceed 
10% of the uncertainty without systematics. 



5.5 Expected Beam Uncertainties 

Before quoting and discussing the allowed levels of differntial ellipticity, gain and beamwidth 
it is instructive to estimate the uncertainty within which these beam parameters will be 
determined from a beam- calibration with a nearly black-body point-source such as Jupiter 
[Tp ^ 200K, 6p ^ 0.5 arcmin). By Wiener filtering a map of the observed source one expects 
to recover the source image with a signal-to-noise level 



Nj J P(l) (27r)2' 

where S{1) is the Fourier transform of the point source and P{1) is the instrumental noise, 
i.e. the {S/N^ is the ratio of the signal and noise power-spectra integrated over all accessible 
multipoles in the experiment in question. The higher S/N the smaller are the uncertainties in 
the recovered beam parameters. The Fourier transform of the convolved point source reads 

T;^« = (1 + g)fpe-'2^''''-'^^>y-'^-P (36) 

where we assume an elliptical gaussian beam with principal axes and ay, gain 1 + g and 
pointing p. This results in 
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V'- (37) 



Since the pointing merely adds a phase to the beam function it drops from the expression 
for S/N. Similarly, S/N is also independent on the beam rotation angle since temperature 
measurements are insensitive to e. t] is an experiment-specific optical-efficiency parameter 
(< 1). Therefore, the following procedure, which is based on S/N considerations, will be used 
to determine the uncertainty of e, /i and g only. To determine these uncertainties we require 
that varying the beam parameters will result in signal changes smaller than the noise 

{S/Nf {S/Nf + 1. (38) 

This condition readily yields the uncertainty in beam parameters 

A(e2) = A(/i2) = A(/) = [S/iNv)]-'. (39) 

For each 'experiment' of the 20 considered here we plug in Af, and 6^ and obtain rjN/ S. Table 
3 summarizes rjN/S in %-units (as in Miller, Shimon & Keating 2008). The color-coded 
threshold values reported in Tables 4, 6 & 7 (differential gain, beamwidth and ellipticity. 
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respectively) should be compared with those of Table 3. Values in green are those which meet 
the requirements from the uncertainties (reported in Table 3) by better than factor 20. Blue 
figures are those which meet the fundamental uncertainties specified in Table 3 in case the 
optical efficiency is 77 = 1 but fail to do so if it is as low as r/ = 0.05. 



5.6 Conclusions 

The upper limits we obtained on the allowed range of beam mismatch parameters for given 
experiments and given arbitrarily-set tolerance levels on the parameter bias and uncertainty, 
constitute conservative limits in the treatment of systematics but on the other hand they 
neglect potential confusion sources with lensing in the worst case scenarios as we explain 
below. It may be the case that a few of the systematics studied here may be fully or partially 
removed. This includes, in particular, the first order pointing error which couples to the 



dipole and octupole moments of non-ideal scanning strategies (see |47|). By removing this 
dipole during data analysis the effect due to the systematic first order pointing error (dipole) 
may drop dramatically. We made no attempt to remove or minimize these effects in this work. 
Our results highlight the need for scan mitigation techniques because the coupling of several 
beam systematics to non-ideal scanning strategies result in systematic errors. This potential 
solution may reduce systematics, which ultimately propagate to parameter estimation, and 
affect mainly the parameters to which the B-mode polarization is sensitive. A brute-force 
strategy to idealize the data could be to remove data points that contribute to higher-than- 
the-monopole moments in the scanning strategy. This would effectively make the scanning 
strategy 'ideal' and alleviate the effect of the a priori most pernicious beam systematics. 
This procedure 'costs' only a minor increase in the instrumental noise (due to throwing 
out a fraction of the data) but may potentially reduce the most pernicious reducible beam 
systematic, i.e. the first order pointing error ('dipole' effect). 

Our results are summarized in Tables 4-8, where we list upper limits on the allowed beam 
systematics (differential gain, pointing, beamwidth, ellipticity and rotation, respectively) for 



various combinations of CMBPOL noise and resolution parameters (in units as in 48|, i.e. 
allowed differential gain, pointing, beamwidth, ellipticity and beam rotation are given in 

\J~^%-, y^arcsec, % and degrees, respectively) based on the requirement that the 

bias induced in the tensor-to-scalar ratio r (first value) and (second value) do not exceed 
the 10% level. Before describing specific results for the various systematics and cosmological 
parameters r and My we comment that all threshold values found from our analysis for the 
differential gain, beamwidth and ellipticity are larger than the beam uncertainties obtained 
in section 5.5 (Table 3) and therefore the beam systematics meet the requirements. 

As for the pointing and beam rotation; these are unconstrained by the beam calibration 
with unpolarized point-source as was illustrated in section 5.5. However, as explained above, 
the effect of pointing may be harnessed by removing non-ideal moments of the scanning 
strategy and addressing the effect of polarization-mixing will require an accurate measurement 
of polarization direction. Tables 4 and 5 show the tolerance for the differential gain and 
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pointing respectively, subject to the bias in r and M^. Our discussion begins with r. As 
expected, when the sensitivity and resolution of the experiment increase - the bounds on the 
allowed systematics are more demanding. Tables 6 and 7 refer to the differential beamwidth 
and ellipticity. Here the dependence on sensitivity is as before but the allowed systematics 
actually increase as the angular resolution improves. The reason for this behavior is simple; 
both differential beamwidth and ellipticity scale as second gradients of temperature. As a 
result, they steeply rise as a function of multipole number and effectively peak on scales 
smaller than the beamsize. This implies, for example, that if we are interested mainly in 
the tensor-to-scalar ratio r we should consider having our beams very narrow so as to push 
the systematic signal to small angular scales, beyond the inflationary peak at ~ 2°. Table 8, 
which describes the allowed beam rotation based on the requirement on r is consistent with 
the general picture we saw with the differential gain and pointing (Tables 4 and 5); increasing 
resolution and sensitivity implies stronger limits on the allowed systematics levels. 

For My forecasts the picture is more complicated mainly due to several competing ef- 
fect and the fact that lensing extraction benefits most from few arcminute scales and few 
of beam systematics peak on sub-beam scales. Higher sensitivity experiments (with 1 and 
2/iK — arcmin) exhibit interesting behavior: the most stringent limits come from 5 '-20' reso- 
lution experiments. This is where the B-mode from lensing peaks and since these low-noise 
experiments are sensitive to the lensing signal even small systematics might potentially bias 
the inferred My. 

Table 5 describes the allowed pointing levels. A low sensitivity experiment (6 /iK-arcmin) 
allows increasing levels of differential pointing as we increase the resolution. However, as 
we increase the sensitivity the most stringent constraints come from experiments with ~ 10' 
resolution as the lower noise level allows to 'see' larger portions of the B-mode lensing signal 
and therefore the allowed systematics are relatively smaller. Now, for a given resolution; 
moderate-low resolution experiments (10'-30'), which do not target the peak of the B-mode 
lensing signal even with systematics- free experiment, allow increasing levels of pointing as we 
'turn on' the instrumental noise. When the resolution is relatively high (3'-5') we still obtain 
increase of allowed systematics with increasing noise but this increase is not monotonic, 
rather - there is a distinguishable 'dip' at around 4/xK-arcmin: a possible explanation is 
that there are two competing effects in action. The first is that increasing the instrumental 
noise naturally allows increasing systematics without significantly affecting the uncertainty 
on inferred parameters. On the other hand increasing the instrumental noise limits lensing 
extraction and therefore relatively increases the weight of the information contained in the 
primary signal. The latter also suffers from contamination of E- and B-modes by beam 
systematics. The interpretation of the ellipticity constraints (Table 7) is similar to that of 
Table 6. Finally, the 'undulations' in the allowed beam rotation, as shown in Table 8, merely 
reflect the relative ratios of primordial E- and lensing-induced B-mode at different multipoles 
accessible at the various resolutions. 

As seen from the above, the allowed beam asymmetry parameters are non-trivial functions 
of both sensitivity and angular resolution as well as the cosmological parameter in question; 
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Table 1: Definitions of tlie parameters associated witli the systematic effects. Subscripts 1 and 2 
refer to the first and second polarized beams of the dual beam polarization assumed in 
this work. 

r or My. They strongly depend on the tolerance criterion (whether it is r that depends 
on the primordial B-mode signal at angular degree scales or My which is extracted largely 
from the higher multipole regime - lensing extraction from the few-arcmin B-mode signal). A 
possibly important factor in this context is the typical angular scales of these various types 
of systematics (i.e their /-dependence). 

Throughout this study we invoked the standard quadratic estimators for lensing recon- 
struction by Hu & Okamoto - this is allowed since we assume the scanning strategy is 
uniform across the sky (and as a result it induces no new typical scale and hence no new 
non-gaussianity) . In practice however, this need not necessarily be the case; there are new 
non-gaussianities induced by the coupling of scanning strategy to the underlying sky and 
beam asymmetry. Consequently, the quadratic estimators will be biased and a more thor- 
ough, Monte-Carlo-Based study, needs to be carried out in order to fully address this issue. 
The effect of differential pointing, in particular, may mimic lensing by 'shifting' features in 
the polarization maps and since its leading order contribution depends on the coupling of 
temperature anisotropy, beam pointing and scanning strategy, the later will cause a mode- 
mode coupling if it is non-uniform, inducing non-gaussianities. However, this very property 
of coupling to scanning strategy can be used, in principle, to remove it, at least partially. In 
addition, the data may be uniformized to some degree by throwing out data points which 
contribute to the non-ideal scanning strategy (see section 5.4). 
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Table 2: Tlie scaling laws for the systematic effects to the power spectra Cf^, Cf^, Cf and Cj^ 
assuming the underlying sky is not polarized (except for the rotation signal where we 
assume the E, and B-mode signals are present) and a general, not necessarily ideal or 
uniform, scanning strategy. The next order contribution ( 10% of the 'pure' temperature 
leakage shown in the table) is contributed by Cf^. It can be easily calculated based on 



the general expressions in 471] where the definitions of z, e, etc., are also found. For 
the pointing error we found that the 'irreducible' contribution to B-mode contamination, 
arising from a second order effect, is extremely small and therefore only the first order 
terms (which vanish in ideal scanning strategy) are shown. The functions /i and /2 are 
experiment-specific and encapsulate the information about the scanning strategy which 
couples to the beam mismatch parameters to generate spurious polarization. In general, 
the functions /i and /2 are spatially-anisotropic but for simplicity, and to obtain a first- 
order approximation, we consider them constants in general. In the case of ideal scanning 
strategy they identically vanish. The exact expressions are given in 47 1. 
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Table 3: Epected uncertainty in the beam parameters e, g and (= r]N/S) from beam calibration 
with point-like source Jupiter {6p ~ 0.5 arcmin), all given in % units, as a function 
of instrument sensitivity and beamwidth (as described in section 5.5). rj (< 1) is the 
experiment-specific optical efficiency; this parameter encapsultes our current ignorance of 
the experiment optics. 
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30' 


0.00183 
0.021 


0.00265 

0.050 
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0.033 


0.00536 
0.049 



Table 4: Tolerance levels for differential gain as a function of instrument sensitivity and beamwidth 
(set by requiring that the fractional error induced in the inferred r (first) and Mi, (second) 
do not exceed 10%). We assume here worst-case-scenario scanning strategy. 
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1.10 
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0.76 
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0.79 
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0.29 


0.09274 
0.62 
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1.60 


0.30476 
1.10 


20' 


0.04853 
0.31 
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0.59 


0.21992 
1.30 
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2.20 


30' 


0.05876 
0.38 


0.12186 
0.70 


0.24634 
1.50 


0.36842 
2.40 



Table 5: Tolerance levels for differential pointing as a function of instrument sensitivity and 
beamwidth (set by requiring that the fractional error induced in the inferred r (first) and 
My (second) do not exceed 10%). We assume here worst-case-scenario scanning strategy. 
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0.015 
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0.032 
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Table 6: Tolerance levels for differential heamwidth as a function of instrument sensitivity and 
beamwidth (set by requiring that the fractional error induced in the inferred r (first) and 
Mjj (second) do not exceed 10%). We assume here worst-case-scenario scanning strategy. 
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0.54 


20' 
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0.02526 
0.087 
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0.17 
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0.25 
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0.00605 
0.029 
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0.054 
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0.10 
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0.14 



Table 7: Tolerance levels for differential ellipticity as a function of instrument sensitivity and 
beamwidth (set by requiring that the fractional error induced in the inferred r (first) 
and Ml, (second) do not exceed 10%). 
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Table 8: Tolerance levels for beam rotation as a function of instrument sensitivity and beamwidth 
(set by requiring that the fractional error induced in the inferred r (first) and Mi, (second) 
do not exceed 10%). 

6 Discussion 

Gravitational lensing imprints the large-scale potentials along the line-of-sight to last scat- 
tering on the observed CMB and generates a guaranteed B-mode signal. The most powerful 
technique for using this extra information is the quadratic estimator formalism, which ex- 
tracts the lensing signal in the form of a noisy map (f)£m of the lensing potential. For the 
sky coverage, noise levels, and beam size expected for CMBpol, this indirect measurement of 
(pim will have high signal-to-noise on a wide range of angular scales, and can be a source of 
cosmo logical information which is complementary to the primary CMB. It is possible to place 
constraints on "late universe" parameters such as neutrino mass (Fig. [5]), the dark energy 
equation of state (Fig. E]), and curvature (Fig. [7]) from the CMB alone. 

In addition to being a source of cosmological information, gravitational lensing is also 
a contaminant of the gravitational wave B-mode signal on large scales. The experimental 
requirement for lensing to be a limiting source of uncertainty is quite stringent: the instru- 
mental sensitivity must be ~ 5 /iK-arcmin or better, and contamination from foregrounds and 
instrumental systematics must also be controlled to better than this level. However, if these 
requirements can be met, CMB experiments will enter a regime in which large-scale B-mode 
measurements are intimately connected with small-scale measurements of the lensing B-mode: 
the only possibility for further improvement in a{T/S) will be to use "delensing" techniques 
which use the B-modes on small scales to reduce the level of lensing contamination. The im- 
provement in ^(T/ S) from delensing will depend on the instrumental noise level and beam. 
Since both instrumental sensitivity and beam width are primary drivers of the total cost 
and complexity of a mission, weighing the tradeoffs is likely to be a complex question when 
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designing experiments, particularly since foregrounds and instrumental systematics will also 
be considerations. Our main result (Fig. ([H])) shows the dependence of the residual B-mode 
noise level on large scales, as a function of the noise and beam, to help in this design decision. 
"External" delensing of the gravity wave B-mode, either via lens reconstruction from small- 
scale CMB temperature anisotropics or large-scale structure, is not a promising approach; we 
establish "no-go" theorems ( §3.2[ §3.3p showing that the improvement in a{T/S) is minimal, 
even under optimistic simplifying assumptions. 

We have done a detailed analysis of the contamination to the lens reconstruction expected 
from polarized foregrounds, and concluded (§4]) that currently favored foreground models do 
not predict that residual foregrounds will be a significant source of bias. Our r and M,^ 
forecast in the presence of beam systematics illustrates that the five types of systematics 
considered here will not significantly bias either the B-mode constraint on (T/S) or the lens 
reconstruction constraint on neutrino mass. We conclude that neither foregrounds nor beam 
systematics are likely to be a limiting factor to the promising science that lies ahead, as 
future generations of experiments probe gravitational lensing through the CMB polarization 
on small angular scales. 
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A Methodology: lens reconstruction and delensing 



A.l Lens reconstruction from CMB temperature and polarization 

We will construct the lens reconstruction estimator (j)£m and its noise power spectrum in a 
uniform way which applies to both temperature (used in §3.21) and polarization. The notation 
in this appendix follows 37 . 

In the presence of a nonzero lensing potential, the CMB two-point function acquires off- 
diagonal (i.e. i ^ m ^ m!) correlations. To lowest order in 0, these take the form 



nil 



(40) 



where X, F G {T, B}. (The two-point function in Eq. fl4Up is the most general form which 
is linear in and satisfies overall rotation invariance.) The F couplings are given by: 
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The estimator is constructed as follows. Assume signal + noise power spectra 
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The minimum variance unbiased estimator and its noise power spectrum are given by: 
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(Here and elsewhere in this appendix, we give expressions in harmonic-space form, but we 
note that for practical evalua tion it is necessary to rewrite them in a computationally efficient 
position-space form, see e.g. |l23l |.) 

This construction assumes that both temperature and polarization are combined into a 
single minimum- variance estimator (pim-, but temperature-only and polarization-only recon- 
structions are the special cases {Nf^, Nf^} and Nf^ — > respectively. 



A. 2 Delensing 

In the gradient approximation, the lensed B-mode is given in terms of the unlensed 13-mode 
and the lens potential by: 

B _ -pEB E* </>* f!-n\ 

where r^ f„^ w as defined in Eq. (1451) . In this notation, the lensed B-mode power spectrum is 
given by 
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^j:Kuc^.'cr (51) 



To perform delensing, we construct an estimator af^ for the lensing B-mode, given the recon- 
struction (pirn described in the previous subsection, and a (noisy) observation of the E-mode 
af^. Heuristically, the estimator is constructed by simply substituting the Wiener-filtered 
Elm and (ptm into Eq. (jSO]): 

fiEE E* r<'t>4>'7* 

-pfB ^ yEB ^li "fimi <P£m (^r)\ 

"'i2m2 2-^ ^ lll2i (JEE _^ J\JEE ^4>4> , -p^clxj) 

iimiim ^1 ^1 £ ' £ 

The delensed B-mode power spectrum, i.e. the power spectrum of the residual field {ciirn~^fm)y 
is given by: 

C,^^ (delensed) = — ^ V \TfU' Cf^Cf" ^ ' I (53) 
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(A more formal derivation of the estimator in Eq. (152|) can be obtained by solving for the 
weights in the estimator which minimize the delensed B-mode power spectrum.) 

The noise power spectrum Nf'^ in Eq. fHSl) and the residual B-mode power spectrum 
(delensed) in Eq. ( 153!) were obtained assuming quadratic lens reconstruction. If we want 



to assume iterative lens reconstruction [33|, l3J], then we simply repeat the calculation of 
iV/"^, C^^( delensed) replacing the lensed B-mode power spectrum (which enters the calcu- 
lation of Nf"^ as a source of noise in Eq. (HSi) ) by the current value of C^-^ (delensed), and 
iterate until C^'^ (delensed) converges. 

We calculate Nf*^, C/'^ (delensed) using the scheme described above when making forecasts 
throughout this report. It is important to note that this scheme is actually an approximate 
forecasting procedure for computing these power spectra and is not exact. For example, a 



53 



complete calculation of the power spectrum C^^(delensed) to lowest order in would include 
an 8-point correlation function containing 6 factors of E and 2 factors of 0. In this language, 
the approximate result in Eq. (IH^ would equal a subset of all the contractions obtained using 
Wick's theorem. To show that our approximate scheme is in fact a good approximation, we 
coinpared forecasts obtained using this scheme with the Monte Carlo simulations of delensing 



m 



15l | for specific choices of noise level and beam, and find good agreement. 



Finally, we discuss forecasts for the parameter uncertainty a{T/S). As mentioned in §31 
the value of cr(T/ S) achievable in a given experiment will depend not only on the instrumental 
noise and beam, but also on sky coverage and loss of modes at low i due to EB mixing 
from survey boundaries, or from projecting out foregrounds. To quantify this latter set of 
complications, we introduce a mode density function {dN^o^cs/di) to represent the number 
of modes that can be measured at a given value of i, and write the parameter uncertainty 
cr(T/ S) with and without delensing as: 



aiT/S) 

no delensing 
^(T/ iS')with delensing 



1 / dNj^odcs 

2^1 d£ 



Cf-^ (tensor) 



Cf ^(lensed) + {a^^J^ 



1 -1/2 



dK 



modes 



di 



Cf^ (tensor) 



Cf ^(delensed) + {a^^J^ 



-1/2 



(54) 
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(We have normalized dN^odcs/di so that dN^odes/ di = {2i +1) for an all-sky survey.) 

Empirically, we find that both the lensed and delensed B-mode power spectra are constant 
at low i to an excellent approximation: 



Cf^ (lensed) 
Cf^ (delensed) 



I '^lensed J 

( B ^ 
V delensed/ 



{i < 100) 
2 {£ < 100) 



(56) 
(57) 



(To quantify this better, the difference between each power spectrum and its best-fit constant 
approximation is below the cosmic variance limit for detectability using only multipoles with 
£ < 100.) From this and Eqs. ([54D, §B it follows that: 



a{T/S) 



no delensing 



'^knsed)^ + ('^inst; 



delensing 

We have used this simplication throughout 



l^felensed)^ + K^^ 



(5^ 



inst / 
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B Methodology: Fisher forecasts 



As discussed in the text, CMB lensing reconstruction allows us to add an additional source 
of information to the usual CMB temperature and polarization fields — namely, the lens 
reconstructed deflection field, d™. The deflection field contains information about late time 
geometry and structure in the universe and helps break the angular diameter distance degen- 
eracy in the CMB. 

For the purpose of Fisher Matrix calculations, it is profitable to assume that the lens 
reconstruction has been used to de-lens the temperature and polarization fields, yielding four 
Gaussian independent variables {T[^ , E]^ , BJ" , d]^} , which are the unlensed CMB fields plus 
the deflection modes. If we further assume, as we have done here, that the fiducial model has 
no primordial B mode, then the de-lensed B-mode is purely noise. Omitting the latter from 
the Fisher calculation, the data covariance matrix reads, 

cr + 

Cf^ + Nf^ I (59) 

where the Cf"^ 's are the unlensed power spectra and the N^^^s denote noise power spectra. 
The deflection field power spectrum N^"^ can be computed in the context of a quadratic 




estimator for the deflection field, as prescribed in |123 



Under these assumptions, the Fisher Matrix can be simply written as, 

(60) 



Fij = ^^^^ /sky Trace 



dpi dpj 



where Pi denotes the i-th cosmological parameter and the lower bound on the error on pi after 
marginalization over all other free parameters is given by. 



For the cases considered in this report we considered a standard 6 parameter ACDM cos- 
mology, parameterized via {VLi,h? ^Vt^yih'^ ^OaiT^Hs^ A^} and extended it to include a massive 
neutrino density as a fraction fy of the total dark matter density, a dark energy equation of 
state w, and a curvature energy density VLk- Note that we chose 6a, the angular scale of the 
sound horizon at recombination, as a parameter, rather than I^a. This is crucial because 9 a 
is the observed quantity and should be kept constant when evaluating the derivatives with 
respect to other parameters. We ensured the convergence of the Fisher calculations by re- 
peating all exercises after halving the step sizes and making sure that the derivatives and 
constraints on parameters remain effectively unchanged. We chose a fiducial model given by 
{VLi,h^, ^^DM/i^ Oa, r, n„ (10^ A,), w, n^} = {0.023, 0.121, 0.010464, 0.11, 0.96, 2.453, 0.008, -1, 
with one massive neutrino species. 
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